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ABSTRACT 



In this letter we present an overview of the rich population of systems with 
multiple candidate transiting planets found in the first four months of Kepler 
data. The census of multiples includes 115 targets that show 2 candidate planets, 
45 with 3, 8 with 4, and 1 each with 5 and 6, for a total of 170 systems with 408 
candidates. When compared to the 827 systems with only one candidate, the 
multiples account for 17 percent of the total number of systems, and a third of 
all the planet candidates. We compare the characteristics of candidates found in 
multiples with those found in singles. False positives due to eclipsing binaries are 
much less common for the multiples, as expected. Singles and multiples are both 
dominated by planets smaller than Neptune; GQlg percent for singles and 86^5 
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percent for multiples. This result, that systems with multiple transiting planets 
are less likely to include a transiting giant planet, suggests that close-in giant 
planets tend to disrupt the orbital inclinations of small planets in flat systems, 
or maybe even to prevent the formation of such systems in the first place. 

Subject headings: planetary systems 



1. INTRODUCTION 



Although it was anticipated that NASA's K epler mission cou l d find systems with more 



than one planet transiting the same host star (IKoch fc Boruckil Il996l : iHolman fc Murray 



20051 ). the rich harvest of candidate multiples that appeared already in the first four months of 



Kepler data, cau ght all of us on the K epler Science Team by surprise. The first announcement 



Steffen et al. 


2010 


) was 


Jsorucki et al.l 


2011a) 



with information that pro vides additional constraints o n the characteristics of the planets 



and even their host stars (IRagozzine fc Holman.l l201ll ). as illustrated by two examples of 



multiple planet systems exhibiting transit time variations that constrain the masses o f the 



planets: Kepler-9 w ith three transiting pl anets (IHolman et al.ll2010l : iTorres et al.ll201ll ). and 



Kepler- 11 with six (ILissauer et al.l 1201 lal ). In this Letter we present an overview of the full 
population of multiples that show transits in the first four months of Kepler data. Note 
that we have not attempted to correct for the probability that planetary orbits are properly 
aligned to show transits, nor for the dependence of transit detectability on various noise 
sources. Instead we have chosen to compare singles with multiples in ways that should 
minimize these biases. 



KEPLER OBJECTS OF INTEREST 



Kepler targets that show features in their light curves that might be due to transits 
are designated "Kepler Objects of Interest" (KOIs). The KOI numbering convention is 
that the digits before the decimal point specify a unique target, and the two digits after 
specify a planet candidate, in the order that it was identified for that target. There is no 
simple description of how KOIs were identified, because the procedures evolved considerably 
as the data improved and the team gained experi e nce. T he general approach used for the 
identification of KOIs is described by lBorucki et al.l (1201 Ibl ). Here we present some additional 
details, with special emphasis on the procedures that were used to identify candidates in 
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multiples. 



Initially, KOIs were identified by visual inspection of light curves for candidates identified 
by the Kepler pipeline using the Transiting Planet Sea rch (TPS; Jenkins et al.l |2010| ) on 



individual quarters of data. The Data Validation (DV; IWu et al. 



2OI0I ) reports from the 



pipeline were then used to identify false positives involving centroid motion during dimmings 
and also to identify additional candidates. This effort resulted in nearly 1000 KOIs and 
somewhat less than 100 systems of multiple candidates. 

The next major release of the Kepler pipeline will stitch quarters together, so that TPS 
and DV can work on light curves from multiple quarters. As a stopgap, a stand-alone tool 
for analyzing multiple quarters was developed by Jason Rowe. Starting with the calibrated 
(raw) time series, sections were excised that showed instrumental artifacts, such as gaps 
due to safe modes of the spacecraft and subsequent thermal settling. The light curves were 
next detrended with a high-pass filter (to reduce sensitivity to instrumental drifts and long- 
term stellar var iability) and then were searched for transits using a version of the Box Least 
Squares (BLS; iKovacs. Zucker. fc Mazehl 120021 ) algorithm. Multi-quarter data for nearly 
180,000 targets were searched for transits (some targets were observed for only one or two 
quarters). Any event that was detected above a 3-a threshold was sent to routines that 
attempted to fit a planetary-transit model, adopting the stell ar parameters (Tes, logs', 
i?^) from the Kepler Input Catalog (KIC; iBrown et al.l 1201 if ). Plots of the successful fits, 
about 25,000 in all, were then inspected visually. This effort led to about 600 additional 
KOIs, with nearly 100 of them in multiples. 

The stopgap multi-quarter pipeline was then run again on the earlier set of KOIs, after 
removing the sections of the light curves affected by the previously identified transits, to look 
for additional candidates. This process was iterated until no more candidates were found. 
This effort identified more than 100 new candidates in multiple systems, in addition to the 
candidates that had been identified previously using the quarter-by- quarter analysis. 



3. FALSE POSITIVES 

KOIs are reviewed from time to time by the Kepler team, to determine which ones 
should be prioritized for additional follow-up observations of various types, and which ones 
are likely false positives that can be retired to the inactive list. In the paper summarizing 
the characteri s tics of t he planet candidates identified in the first four months of Kepler data. 



Borucki et al.l ( 1201 Ibl ) present a list in their Table 4 of 498 KOIs that had been identified 



as false positives and were no longer considered to be viable planet candidates. In Table 1 
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we summarize the number of false positives compared to the number of surviving candidates 
among the KOIs, with a separate accounting for the singles and multiples. 

More than half of the false positives (59%) resulted from an "active pixel offset" (APO). 
This test uses a difference image analysis to show that during transit-like events the image 
is significantly displaced from the target position and is star-like, indicating contamination 
by a faint background eclipsing binary or by the wings of the PSF from a nearby bright star 
encroaching on the edge of the target aperture. Most of the remaining false positives also 
involved eclipsing binaries, and were identified either by features in the Kepler light curves, 
such as secondary eclipses or ellipsoidal variations or eclipse-time variations (34%), or by 
large variations observed in the radial velocities (5%). Eleven of the early KOIs were judged 
to be photometric false alarms, based on additional data from subsequent quarters. 

The difference in the rate of false positives for singles compared to multiples is striking; 
for the singles the rate is 37 percent (486 false positives compared to 827 survivors), but for 
the multiples the rate is only 3 percent (12 false positives in 6 systems, compared to 408 
surviving planets in 164 systems). This large difference is expected, because the APOs for 
the singles are the result of chance alignments of eclipsing binaries with the full target list of 
nominally 150,000 stars, while the APOs for the multiples come from chance alignments with 
systems that already show transit-like events. Lumping together all the false positives among 
singles due to eclipsing binaries gives a rate of (288+164+23)/150,000 = 0.0032. Assuming 
that the KOIs have been drawn from the same parent population as all 150,000 targets, 
the expected number of doubles involving a planet and a false positive due to an eclipsing 
binary is roughly 827 x 0.0032 = 2.6, while the number involving two eclipsing binaries and 
no planets is roughly 486 x 0.0032 = 1.5. The key assumptions here are that the probability 
that a target image is contaminated by an accidental alignment with a background eclipsing 
binary is the same, whether the target already shows a transit-like event or not, and that 
the false-positive probability for a KOI is independent of the depth of the transit-like event 
in its light curve (multiples mostly show shallower dips). By the same line of argument, the 
number involving two planets and an eclipsing binary is only 115 x 0.0032 = 0.4. Higher 
order coincidences are correspondingly less likely. The probability that a double consists 
of an accidental alignment of two unrelated singles seems less likely than the two eclipsing 
binary case, because eclipsi ng binaries are m ore common than singles in the A'e]5/er sample 



by a factor of almost three (jPrsa et al.ll201ll ). These numbers are summarized in Table 1. 



These rough estimates of the expected rates of false positives among the multiples are 
preliminary, because the vetting effort is still unfinished. Furthermore, some types of false 
postives, such as hierarchical triples, are extremely difficult to identify, especially for shallow 
events. Thus we expect that there are still false positives lurking among both the singles 
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and the multiples. Nevertheless, in general it must be true that false positives due to chance 
alignments must be much less common among multiples than singles, and even for singles 
the rate of residual false positive s among the vetted candidates may be as low as 5 or 10 
percent flMorton fc Johnsonll201ll ). 



4. PLANET RADIUS VS ORBITAL PERIOD 

The process of fitting transit models to KOI light curves delivers two primary observable 
characteristics of the candidate planet: the planetary radius, Rp (where we have adopted 
the KIC value for the stellar radius), and the orbital period, P. The plot of these two 
quantities against each other is shown in Figure 1, where the active KOIs in multiples are 
blue and the singles are red. Single planets come in all sizes, but there are relatively few 
giant planets in transiting multiples. The pile-up of giant planets near 3 days is obvious, 
with no corresponding pile-up of planets, either large or small, among the multiples. The 
detection limit is especially clear in the lower right corner of the upper panel, where period 
is logarithmic. The edge of the distribution of detected candidates has a slope of 1/3 as 
expected (the total number of data points during transits varies nominally as P~^/^). 

The distributions of planet radius versus period are shown more quantitatively by the 
histograms in Figure 2, where we have collapsed Figure 1 onto the -Rp and logP axes in 
the upper and lower panels, respectively. The vertical scales for the singles and multiples 
have been normalized so that they both have the same area under their histograms. Planets 
smaller than Neptune dominate both samples, but more so for the multiples; 69^3 percent 
for the singles and 86I5 percent for the multiples. The error estimates for these percentages 
only consider Poisson noise and do not include any contribution from uncertainties in Rp. 
The difference in the radius distributions between the singles and the multiples is highly 
significant; the K-S test gives a probability of 2 x 10^^*^ that they are drawn from the same 
parent distribution. 

The period distributions for singles and multiples are quite similar. To the eye there 
may appear to be a slight shift to shorter periods for the singles, but the significance of this 
difference is not supported by the K-S test, which reports a probability of 10 percent that 
such differences could occur by chance. 

Figure 3 compares the number of singles versus the number of systems that are multiples, 
as a function of effective temperature. Because nearly all of the host stars are on or near 
the main sequence, effective temperature is a reasonable proxy for host-star mass. The K-S 
test reports that the difference between the two distributions is marginally significant, with 
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a probability of 0.008 and D value of 0.1. It appears that singles may be more common 
than multiples around the hotter, more massive stars, while the multiples are more common 
than singles around the cooler, less massive stars. This might be rel ated to the tendenc y 
for close-in giant planets to be less common around low-mass stars (j Johnson et al.ll2010l ). 
and/or to the tendency for small planets in the Rp range 2 t o 4 Earth radii to be much more 
common around cool stars in the Teg range 3600 to 4100 K (IHoward et al.ll201ll ). Note that 
this is only a comparison of singles to multiples, and no corrections have been made, either 
for the relative number of targets as a function of T^s or for the probability of detection. 



5. DISCUSSION 



The fraction of single planet candidates that are smaller than Neptune is 69lil percent 
(569/827). The fraction of multiple systems with no planets larger than Neptune is 78^1^7 
percent (133/170). Thus, systems with multiple transiting planets are less likely to include 
a transiting giant planet. If the comparison is restricted to short period planets (P < 10 
days), the difference is particularly striking: the fraction of short-period single candidates 
that are smaller than Neptune is 691^4 percent (279/405), while the fraction of multiple 
systems that contain at least one short-period planet (117 systems) but no short-period 
planets larger than Neptune is 96^9 percent (112/117). One possible interpretation is that 
a close-in giant planet can stir up the orbits of other inner planets in its system, while a 
system of small planets is more likely to preserve the flatness of the disk from which it formed. 
This picture is supported by determinations of the spin/orbit alignment for transiting giant 
planets using the Rossiter-McLaughlin effect. The orbits of some giants are well aligned with 
the rotation of their ho st star, while others show significant orbital inclinations, including 
even retrograde orbits (jWinn et al.ll2010l : iTriaud et al.l 120101 ). Thus there is good evidence 
that some systems have been disrupted from their presumably flat initial configuration. On 
the other hand, there is good evidence that close-in giant planets do not always disrupt 
or prevent the formation of systems with multiple planets. Radial-velocity surveys show 
that ab out 25% of the giant planets are acco mpanied by companions with smaller minimum 
masses (jWright et al.ll2009l : ISchneiden l201ll ). and the actual fraction may be much higher 
due to the radial-velocity detection limit for small planets. 

We observe the rate of false positives due to eclipsing binaries to be much smaller for 
multiples than singles. This is expected, because the number of candidates that show a 
candidate planet or false positive is much smaller than the full list of approximately 150,000 
targets. Thus, the probability th at an eclipsing binary contaminates the light of a multiple 
is much smaller than for a single. iLissauer et al.l (l2011bl ) present some independent evidence 
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that many of the multiples must be systems of planets, in particular the common occurrence 
of periods near mean motion resonance. This reinforces the impression that KOIs in multiples 
are very likely to be planets. 

Why hasn't CoRoT announced any multiples yet? Kepler's better photometric precision 
and longer time series both contribute to the detection of smaller planets with longer periods, 
as is needed to discover fiat systems. Actually, CoRoT may have come very close to detecting 
a multiple transiting system, namely CoRoT-7. The transiting planet in this system, CoRoT- 
7b, is the smallest discovered so far by CoRoT (see Figure 1), but the orbit has a rather 
extreme impact paramete r. Thus additiona l planets in the system, such as the proposed 



second planet CoRoT-7c (iQueloz et al.ll2009l ). would be less likely to transit also 



Transit time variations for planets in multiple systems promise to be an important tool 
for constraining the masses of planets that are too small to be detected with current radial- 
velocity techniques. These constraints improve with longer time series, which is a good 
argument for extending the Kepler mission and continuing to monitor the most promising 
multiple systems. This approach may b e able to confirm rocky planets in the Habitable 



Zones of Kepler targets (IFord et al.ll201ll ). 



We thank the entire Kepler team for all the hard work that has made these results possi- 
ble. Funding for this Discovery Mission is provided by NASA's Science Mission Directorate. 
We give special thanks to the anonymous referee for insightful and timely feedback. 
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Table 1. Census of False Positives 



Number 


Description 


Fraction 


1733 


KOIs identified for 1489 targets 




1235 


Survivor KOIs for 997 targets 


1235/1733 = 0.713 


827 


Survivors in singles 


827/1235 = 0.670 


408 


Survivors in multiples 


408/1235 = 0.330 




170 multiple systems 


170/997= 0.171 




115 systems with 2 KOIs 


230/408 = 0.563 




45 systems with 3 KOIs 


135/408 = 0.330 




8 systems with 4 KOIs 


32/408 = 0.078 




1 system with 5 KOIs 


1/408 = 0.002 




1 system with 6 KOIs 


1/408 = 0.002 




raibe pooiLiveo 


^UOj 1 / oo — U.Zoo 


486 


False positives in singles 


486/498 = 0.976 




Rate per KOIs in singles 


486/(486 + 827) = 0.370 


288 


Active pixel offsets (APO) 


288/486 = 0.592 


164 


Light curve evidence 


164/486 = 0.337 


23 


Spectroscopic binaries 


23/486 = 0.047 


11 


Photometric false alarms 


11/486 = 0.023 


12 


False positives in 6 systems 


12/498 = 0.024 




Rate per KOIs in multiples 


12/420 = 0.029 


5 


APOs in 3 systems 


3/170 = 0.017 


5 


Light curve evidence 


5/420 = 0.012 


2 


In one spectroscopic binary 


1/170 = 0.006 



10 

Period [days] 



100 
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Fig. 1. — Planet radius versus orbital period. The bottom panel shows periods shorter 
than 60 days using a linear scale, the top panel uses a log scale to show the entire range of 
periods. Planet candidates in singles are plotted in red, and in blue for those in multiples. 
Kepler finds very few giant planets in systems of multiple p lanets. This c o nclusio n is not 
affected by the rather large upper limit that was adopted by iBorucki et al.l ( 1201 Ibl ) for Rp. 
The published CoRoT planets are plotted in green. Only one of these, CoRoT-7b, is smaller 
than Neptur ie, and the radial- v elocity observations suggest that it may have a non-transiting 
companion (IQueloz et al.l 120091 ) 
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Fig. 2. — Histograms for the number of planet candidates versus planetary radius and period. 
Singles are shown in shaded red, multiples in cross-hatched blue. The vertical scales for the 
singles and multiples have been normalized so that they both have the same area under 
their histogram. Planets smaller than Neptune dominate both samples, but more so for the 
multiples; 69% for the singles and 86% for the multiples. 
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Fig. 3. — Number of systems versus effective temperature of the host star, which is a proxy 
for the stellar mass. Single planets are relatively more common around the hotter stars, 
multiple planets are more common around the cooler stars. No corrections have been made 
for the relative number of targets as a function of Tefj , or for the probability of detection. 



